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INTRODUCTION 


Decompression  sickness  (DCS)  In  divers  Is  presumed  to  result  somehow  from 
the  accumulation  of  dissolved  Inert  gases  In  tissue,  which  must  occur  during 
periods  of  Increased  ambient  pressure,  l.e.  during  dives.  It  seems  likely 
that  bubbles  produced  during  or  after  decompression  constitute  at  least  one  of 
the  direct  causes  of  DCS  symptoms.  The  risk  of  decompression  sickness, 
therefore,  may  depend  In  part  on  the  solubilities  of  Inhaled  gases  In 
biological  fluids  and  tissues.  Quantitative  prediction  of  DCS  risk  may  be 
aided  by  the  availability  of  measured  values  of  these  solubilities  over  a 
range  of  partial  pressures  relevant  to  divers. 

The  techniques  employed  for  measuring  llquld/gas  equilibria  usually 
Involve  removing  an  aliquot  of  one  of  the  phases  after  equilibration.  The 
composition  of  this  aliquot  Is  assayed  and  the  composition  of  the  other  phase 
Is  then  computed  by  mass  balance.  The  necessary  equilibration  time  Is 
determined  empirically.  Most  of  the  various  experimental  methods  have 
differed  fundamentally  only  In  the  type  of  assay  chosen:  gas  chromatography 
(1-8),  mass  spectrometry  (9-11),  radioassay  (12,13),  and  dllatometry  (14,15). 
The  dllatometrlc  assay  consists  of  reducing  the  pressure  on  a  liquid  aliquot 
and  then  holding  It  at  that  reduced  pressure;  this  Is  done  by  expanding  the 
vessel  containing  the  aliquot  as  gas  desorbs  from  the  liquid  phase .  The 
expansion  can  be  accomplished  by  adjusting  the  piston  In  a  syringe  barrel  that 
constitutes  part  of  the  volume  of  the  gas  phase  above  the  liquid  aliquot.  One 
computes  the  amount  of  desorbed  gas  from  the  necessary  volume  Increase  and 
then  computes  the  composition  of  the  two  phases  by  mass  balance.  This  method 
therefore  requires  two  equilibrations  rather  than  one,  which  provides 
additional  sources  of  uncertainty  in  the  calculation  of  solubility. 

A  simpler  technique  for  measuring  llquld/gas  equilibria  dispenses  with 
the  need  for  separating  an  aliquot  of  material  and  thereby  has  the  potential 
for  greater  precision.  In  this  method  a  liquid  sample  In  a  closed  vessel  Is 
compressed  to  a  known  pressure,  either  by  forcing  gas  Into  the  vessel  or  by 
decreasing  the  vessel's  volume.  After  equilibration,  the  pressure  in  the  5 
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vessel  is  measured  again.  The  amount  of  gas  adsorbed  by  the  liquid  can  be 
calculated  from  the  pressure  change  that  occurs  during  equilibration  (16,17). 

The  published  reports  of  gas  solubility  that  we  have  reviewed  are 
universally  marred  by  the  authors'  failure  to  evaluate  either  the  accuracy  or 
the  precision  of  their  measurements  using  a  formal  error  analysis.  In 
particular,  there  is  no  evidence  that  the  authors  measured  gas  leakage  from 
their  systems  or  that  they  made  a  quantitative  assessment  of  leakage  as  a 
source  of  biased  error. 

We  are  instituting  an  investigation  designed  to  measure  the  solubilities 
of  nitrogen,  helium,  and  other  diving  gases  in  water,  blood,  and  various  solid 
tissues  at  37*0  and  pressures  between  1  and  10  atm.  Solubilities  will  be 
measured  using  the  last  technique  mentioned  above:  the  pressure  change  that 
occurs  in  a  closed  system  during  interphase  equilibration  will  be  measured 
after  an  abrupt  change  in  system  pressure  has  been  effected.  The  experimental 
plan  is  described  in  the  next  section;  it  is  briefly  summarized  as  follows. 

Two  identical  vessels  Immersed  in  a  water  bath  at  37*0  are  bridged  by  a  high- 
precision  differential  pressure  transducer  and  also  by  a  tube  that  bypasses 
the  transducer.  One  vessel  contains  a  measured  amount  of  water  or  tissue 
along  with  gas,  and  the  other  contains  only  gas.  Prior  to  an  experiment,  with 
the  tube  between  the  vessels  opened,  sufficient  time  is  allowed  for  1)  the 
contents  of  both  vessels  to  reach  37“C,  2)  the  gas  phase  of  both  vessels  to 
reach  the  dew  point  through  gradual  evaporation  of  some  of  the  water,  and  3) 
the  condensed  phase  to  equilibrate  with  the  gas  phase  at  37*C  and  known 
initial  pressure  (equal  to  the  ambient  pressure  if  the  system  is  in 
communication  with  air;  otherwise,  equal  to  the  gauge  pressure  plus  ambient 
pressure).  The  initial  conditions  of  the  experiment  are  thereby  established. 
At  the  start  of  the  experiment  the  two  vessels  are  simultaneously  pressurized 
and  the  vessels  then  are  closed  off  from  each  other.  Sufficient  time  for 
reequilibration  is  allowed,  and  the  pressure  difference  between  the  vessels  is 
read  from  the  differential  transducer.  The  amount  of  gas  dissolved  in  the 
liquid  or  tissue  is  calculated  from  this  pressure  differential. 
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In  this  report  the  experimental  method  and  the  means  of  analyzing  the 
data  are  discussed.  Since  chemists,  physicists,  and  engineers  commonly 
consider  gas/liquid  equilibria  in  terms  of  Henry's  Law  or  deviations  from 
Henry's  Law,  we  include  a  detailed  derivation  for  the  equation  to  calculate 
the  Henry's  Law  coefficient  from  the  results  of  the  experiment.  We  also  show 
how  to  compute  the  Ostwald  coefficient,  since  it  has  broader  application  in 
biological  systems  than  does  the  Henry's  coefficient.  Sources  of  error  and 
means  of  limiting  them  are  discussed  quantitatively.  This  intensive  error 
analysis  is  necessary  in  order  to  ensure  the  accuracy  of  the  planned 
solubility  measurements . 
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Experimental  method 

The  experimental  apparatus  Is  diagrammed  in  Fig.  1.  Two  vessels  and  a 
system  of  stainless  steel  tubing  are  immersed  in  an  insulated  water  bath 
maintained  at  37  ±  0.05'C  by  exchange  with  a  heated  supply  bath.  The 
temperature  of  the  water  bath  is  monitored  by  a  thermistor  whose  calibration 
is  updated  periodically  using  a  mercury  thermometer  certified  by  the  National 
Bureau  of  Standards.  The  immersed  portion  of  the  system  is  sjmmetrical  about 
a  differential  pressure  transducer  (T1  in  the  diagram)  and  a  valve  (VI  in  the 
diagram)  through  which  gas  can  be  fed  from  a  cylinder.  Valves  V2  and  V3  allow 
communication  between  the  vessels  and  valve  VI;  valves  V4  and  V5  allow 
communication  between  the  vessels  and  the  atmosphere.  Each  vessel  is  always 
in  communication  with  the  differential  transducer  and  each  vessel  is 
penetrated  by  a  pressure  transducer  referenced  to  the  atmosphere  (T2  and  T3) . 

All  of  the  pressure  transducers  are  calibrated  using  a  dead-weight 
standard  from  CEC  Instrument  Company  (Pasadena,  California) .  Appendix  A 
contains  a  listing  of  the  components  used  in  this  experiment  and  their 
specifications.  The  procedure  we  will  use  for  measuring  gas  solubility  in  a 
condensed  phase  is  as  follows. 

A  sample  of  known  volume  of  the  condensed  phase  is  loaded  into  what  we 
shall  call  the  "experimental"  or  "X"  vessel,  whose  volume  has  been  measured  by 
the  procedure  discussed  in  Appendix  B.  Flask  "R"  is  a  "reference"  vessel 
containing  no  condensed  phase.  The  condensed  phase  might  consist  of  water, 
blood,  or  an  Isotonic  saline  solution  containing  bits  of  tissue.  In  some 
experiments  the  condensed  phase  will  be  allowed  to  equilibrate  with  air  at 
37*C  and  ambient  pressure  by  leaving  valves  V2,  V3,  and  V5  opened,  and  valves 
VI  and  V4  closed.  If  the  condensed  phase  must  be  equilibrated  with  some  gas 
other  than  air,  the  system  will  be  purged  with  the  desired  gas.  Valve  V5  then 
will  be  closed,  and  the  condensed  phase  will  be  allowed  to  equilibrate  with 
the  gas  phase.  In  either  case,  the  system  is  left  undisturbed  overnight  to 
allow  the  gas  phase  to  become  humidified  at  37*C.  This  time  is  more  than 
sufficient  for  the  gas  in  both  flasks  and  in  the  intervening  tubing  to  closely 
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FIGURE  1 


THE  SOLUBILITY  APPARATUS 


approach  the  dew  point  by  passive  diffusion  of  the  water  vapor  evolved  at  the 
liquid  surface.  This  statement  is  supported  by  drawing  the  analogy  between 
mass  transfer  and  heat  transfer  and  using  published  solutions  for  the 
unsteady- state  heating  of  a  homogeneous  bar.  The  solution  for  an  infinitely 
long  bar  heated  at  one  end  is  especially  useful  because  it  is  certain  to 
underestimate  the  humidification  rate  in  our  finite  system,  and  therefore  we 
can  say  confidently  that  it  provides  a  conservative  estimate  of  the  necessary 
humidification  time.  For  the  details  of  these  solutions  the  reader  is 
referred  to  chapter  11  of  Bird,  Stewart,  and  Lightfoot  (18). 

A  magnetic  stirring  bar  is  used  to  agitate  the  contents  of  the  X  vessel 
during  the  equilibration  period,  in  order  to  minimize  transport  resistance  in 
the  condensed  phase.  After  being  left  overnight  the  system  has  the  following 
conditions;  1)  the  system  is  at  37°C,  2)  the  gas  phase  in  both  vessels  is 
saturated  with  water  vapor,  and  3)  the  condensed  phase  is  equilibrated  with 
the  gas  phase  at  a  known  initial  pressure  (i.e.,  the  barometric  pressure  plus 
the  gauge  pressure) .  These  constitute  the  initial  conditions  for  the 
experiment.  The  barometric  pressure  is  measured  using  a  mercury  manometer 
referenced  to  vacuum.  To  begin  the  experiment,  a  pure,  dry  gas  is  injected 
quickly  into  the  system  by  opening  valve  VI  (valve  V5  is  closed  first  if  it 
had  been  open  overnight).  The  injection  must  be  completed  within  a  few 
seconds  to  ensure  negligible  diffusion  of  gas  into  the  condensed  phase  during 
initial  pressurization.  The  symmetry  of  the  system  should  ensure  that  the  two 
vessels  are  pressurized  equally.  The  tubing  has  a  high  enough 
length- to -diameter  ratio  to  preclude  appreciable  backmixing,  so  that  no 
material  should  leave  the  system  during  the  pressurization.  After 
pressurization,  valves  V2  and  V3  are  closed  simultaneously  to  isolate  the 
flasks  from  one  another.  The  contents  of  the  X  vessel  are  stirred 
continuously  during  the  experiment.  The  system  is  allowed  sufficient  time  to 
reequilibrate  and  then  the  following  information  is  recorded;  (1)  the 
indicated  pressure  difference  between  vessels  R  and  X  and  (2)  the  gauge 
pressure  in  vessel  R.  The  barometric  pressure  is  measured  again  so  that  the 
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gauge  pressure  can  be  converted  to  absolute  pressure.  Next,  valves  V4  and  V5 
are  opened  simultaneously  so  that  both  vessels  are  vented  to  the  atmosphere 
and  the  pressure  across  any  one  of  the  transducers  becomes  zero.  The  offset 
voltage  of  each  transducer  (i.e.,  the  indicated  voltage  at  zero  differential 
pressure)  is  recorded  and  its  value  is  used  to  adjust  the  transducer  reading 
taken  a  moment  earlier,  thereby  minimizing  response  drift  as  an  error  source. 

The  pressure  difference  between  the  X  and  R  sides  is  a  function  of  the 
amount  of  gas  that  dissolves  during  the  experiment.  The  placement  of  the 
differential  transducer  between  two  pressurized  flasks  means  that  the  pressure 
across  it  is  much  lower  than  the  absolute  pressure  on  either  side.  This 
allows  us  to  use  a  differential  transducer  rated  for  a  maximum  differential 
pressure  of  only  1000  mmHg  (1.3  atm)  at  absolute  pressures  of  up  to  10  atm. 

The  injection  of  a  pure,  dry  gas  at  the  start  of  the  experiment  effects  a 
change  in  the  partial  pressure  of  the  injected  gas  only.  Other  gases,  e.g. 
oxygen  and  water,  do  not  undergo  a  change  in  partial  pressure.  In  an  ideal 
system,  where  by  definition  the  solubilities  depend  only  on  the  partial 
pressures  of  the  solutes,  oxygen,  water,  and  any  other  miscellaneous  gases 
will  undergo  zero  net  interphase  transport  during  the  experiment  and  will  have 
no  effect  on  the  final  differential  pressure. 

Part  II  of  'Data  analysis"  will  explain  how  the  information  obtained  in 
the  experiment  described  above  can  be  used  to  calculate  gas  solubility. 
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Data  analysis 


I.  Theory:  Equilibria  of  sparingly  soluble  gases 
Discussions  in  this  section  are  adapted  from  Prausnitz  (19). 

The  analysis  of  phase  equilibria  is  greatly  simplified  by  reducing  the 
pressure -volume -temperature  information  to  a  single  parameter,  the  fugacity, 
which  has  units  of  pressure.  The  fugacity  is  defined  such  that,  at 
equilibrium,  its  value  for  each  chemical  species  is  the  same  in  all  phases. 

Thus,  for  a  two-phase,  two - component  system. 


c  B 

-2 


[la] 

[lb] 


where  superscripts  A  and  B  refer  to  the  two  phases  and  subscripts  1  and  2 
refer  to  the  two  species.  For  solutions  of  gases  in  liquids  equation  [1]  is 
usually  rewritten  in  the  following  form  for  the  solute,  which  we  denote  by  the 
subscript  "2": 

h  ~  *2  Yz  ^  ~  Y2  *2  ^2®  [2] 

where  “  fugacity  coefficient  for  the  solute  in  the  gas  pht^se 
(dimensionless) ; 

yj  “  mole  fraction  of  the  solute  in  the  gas  phase  (dimensionless) ; 

P  -  total  pressure  (mmHg) ; 

Y2  “  activity  coefficient  of  the  solute  in  the  liquid  phase 
(dimensionless) ; 

X2  ~  mole  fraction  of  the  solute  in  the  liquid  phase  (dimensionless) ; 
f2°  -  fugacity  of  the  solute  at  some  arbitrary  reference  state  (mmHg) . 
The  value  of  the  activity  coefficient  y2  depends  on  the  value  of  f2°  at 
the  chosen  reference  state.  For  an  ideal  gas  the  fugacity  coefficient  is 
identically  1.0,  and  for  a  liquid  solution  that  is  ideal  in  the  Henry's  Law 
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sense  Y2  is  independent  of  Xj.  Then,  for  an  ideal  system  equation  [2]  becomes 
"Henry ' s  Law" : 

yj  P  -  H  X2  [3a] 

with  H  -  Y2  ^2°  [3b] 

where  H  is  the  Henry's  Law  coefficient  (mmHg/mol  fraction)  and  is  independent 
of  pressure.  The  Henry's  coefficient  is  empirically  determined.  For  anv 
liquid/gas  system  there  is  a  pressure  below  which  Henry's  Law  is  valid.  On  a 
molecular  level,  at  low  mole  fractions  a  solute  molecule  "sees”  only  solvent 
molecules  and  therefore  its  environment  does  not  change  appreciably  with 
changes  in  the  mole  fraction;  this  accounts  for  the  constancy  of  the  activity 
coefficient  and  the  wide  validity  of  Henry's  Law  for  sparingly  soluble  gases. 

At  sufficiently  high  pressures  we  can  expect  non- ideal  behavior  in  either 
phase  (or  both).  In  the  gas  phase  the  fugacity  coefficient  of  a  component  in 
a  mixture  is  given  by  the  following  exact  relation,  which  is  of  general 
validity: 

RT  In  (fi)  -  /v"l  <  T.v.nj  -  <RT/V  )  ]  dV  -  RT  In  (z)  [4] 

where  R  -  gas  law  constant  (mmHg*  cm^)/(gmol' K) ; 

T  -  absolute  temperature  (*K); 

P  -  total  absolute  pressure  (mmHg) ; 

V  -  volume  of  gas  phase  (cm®)  ; 

n^  -  number  of  moles  of  the  species  of  interest  in  the  gas  phase  (gmol) ; 
nj  -  number  of  moles  of  any  other  species  in  the  gas  phase  (gmol) ; 
z  -  compressibility  factor  -  PV/nRT  (dimensionless). 

To  evaluate  the  integral  in  equation  [4]  we  need  either  (1)  a  state 
function  for  the  solute  (i.e.,  a  P-V-T  equation)  that  is  explicit  in  pressure 
or  (2)  compressibility  data  at  the  temperature  of  interest  that  are  sufficient 
to  permit  numerical  integration.  For  evaluating  the  pressure  dependence  of 
the  solute  fugacity  in  the  liquid  phase,  the  Krichevsky-Kasarnovsky  equation 
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is  valid  if  X2  is  sufficiently  small  and  the  solution  is  well  below  its 
critical  temperature: 


ln(Y2f2°)  -  InCfj/Xa)  -  In  (  H  I  Pj-Pi*)  +  [V2-  (P  -  Pi')]/RT  [5] 


where  (  H  I  P2  -  Pi“  )  -  Henry's  coefficient  evaluated  at  the  vapor  pressure  of 

the  solvent  (misHg/mole  fraction) ; 

V2"  -  molar  volume  of  the  solute  in  the  liquid  phase  at  infinite 
dilution  (cm®/gmol) ; 

P^*  -  vapor  pressure  of  the  solvent  at  the  solution  temperature 
(mmHg) . 

Equation  [5]  has  been  shown  to  agree  excellently  with  the  data  on 
nitrogen  and  hydrogen  in  water  at  pressures  as  high  as  1000  atm.  Combining 
equations  [2],  [4],  and  [5]  we  find  that  we  can  relate  y2P  to  X2  using  the 
following  proportionality  factor: 


y2p/x2  -  - 


exp 


ln(H  P2-Pi*l  +  (vz"  (P  -  Pi*)l/RT 


exp 


(  1/RT  )  /v  (  (  dP/dni)  j  y.nj  "  <  i  dV  -  In  z) 


where  we  define  as  an  "effective"  Henry's  coefficient.  For  ideal  systems 
the  complicated  expression  on  the  right  hand  side  becomes  a  constant  equal  to 
H,  and  the  above  relation  decomposes  to  equation  [3],  Henry's  Law.  Therefore, 
a  plot  of  y2P  versus  X2  is  linear  within  the  Henry's  Law  region  and  curved  at 
higher  pressures. 
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II.  Computation  of  the  Henry's  Law  coefficient 

The  nomenclature  for  this  section  is  listed  in  Appendix  C. 

Initial  conditions  on  the  X  vessel  (the  one  containing  two  phases'): 

Conservation  of  the  inert  gas  at  the  start  of  the  experiment  leads  to 

“  V'  +  "0°"  [6] 

That  is,  just  after  the  pressurization  at  time  -  0,  the  liquid  phase  is  still 
equilibrated  with  gas  at  the  old  pressure  but  the  gas  phase  is  at  the  new 
pressure  Assuming  the  mole  fraction  of  the  solute  is  small  enough  that 

(ni+n2)-ni,  we  can  rewrite  equation  [6]  in  terms  of  observable  quantities  as 

-  [(Vi.  Pl)/(H°-  M^)  ]  (P°-  -  Pi'>*  - 

+  [Vs/(2/^T®*)  ]  -  Pi°^  -  P„th°*)  [7] 

The  Henry's  coefficient  H  can  be  what  we  shall  call  an  "effective  Henry's 
coefficient"  if  Henry's  Law  is  not  actually  valid  at  pressure  P®*.  We  define 
an  effective  Henry's  coefficient  as 

-  Vz  PA2  [8] 

without  making  the  fundamental  assumption  of  Henry's  Law,  which  is  that  H  is 
independent  of  pressure. 

Ve  can  simplify  equation  [7]  somewhat  by  recognizing  that  the  system 
temperature  should  equal  the  water  bath  temperature  both  before  pressurization 
and  a  long  time  after  pressurization,  so  that 

r  -  T°'  [9] 

Therefore,  we  can  discard  the  superscripts  on  T  that  denote  times  0-  and 
•B.  Further,  we  assume  that  gases  other  than  the  test  gas  and  the  solvent 
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vapor  are  present  at  low  enough  partial  pressures  that  they  behave  ideally. 
The  following  then  holds  true  for  their  portion  of  the  gas  phase: 

Poth°^  -  (T°VT°-)  Poth“'  -  (T°Vt-)  Poth”'  [lOa] 

Also,  provided  the  system  pressurization  is  accomplished  quickly  enough 
that  there  is  only  negligible  evaporation  of  solvent  or  condensation  of 
solvent  vapor  between  time  -  0-  and  time  -  0+,  then  we  can  write  an  analogous 
equation  for  the  partial  pressure  of  solvent  vapor: 

P^o+  .  (T°VT°‘)  Pi°'  -  (T°''  /T")  Pi°'  [10b] 

Substituting  equations  [9]  and  [10]  into  equation  [7]  yields 

ntot.i“  -  [(Vl  Pi,)/(H°-  Mi)](P‘>-  -  -  P„th°‘) 

+  [Vg  /(z/^T°^)]  -  [Vg  /(z‘’^T“)](Pi°-  +  [11] 

Final  conditions  on  the  X  vessel:  To  describe  the  contents  of  vessel  X  after 
the  two  phases  have  been  in  contact  for  a  long  time  we  have  the  following 
expressions,  which  are  analogous  to  equations  [6]  and  [7]: 

Htofi'  -  +  V  [12] 

and  nt„t.i“  -  l(\  Pl)/(H-  M^)  ]  (Px"  -  Pi*  -  P^tb*) 

+  [Vg  /(Zx-RT-)](Px-  -  Pi-  -  P„th-)  [13] 

To  allow  for  leakage  from  the  "X"  vessel  we  shall  define  parameter  fx  as 
follows : 


“  n^otmi  /ntot«i 


[14] 


12 


As  defined,  fx  -  1  for  a  perfectly-sealed  vessel  and  fx  <  1  if  there  is  a 
finite  leak  when  the  vessel's  contents  are  compressed  above  ambient  pressure. 
Ue  can  rearrange  equation  [14]  to  obtain  the  following  mass  balance  around  the 
"X"  vessel,  which  incorporates  the  possibility  of  leakage: 

To  the  extent  that  Henry's  Law  correctly  describes  the  behavior  of  gases  other 
than  the  test  gas  and  the  solvent  vapor,  and  leakage  is  absent,  their  partial 
pressures  are  the  same  at  time  -  o»  as  it  at  time  -  0-.  This  is  because  1) 
there  is  no  transport  of  these  miscellaneous  gases  across  system  boundaries 
between  those  times  and  2)  the  system  temperature  is  the  same  at  both  times. 

If  leakage  does  occur,  Henry's  Law  indicates  that  the  proportionality  between 
and  remains  unchanged,  from  which  it  follows  that 

Poth“  -  fx  Poth*’’  [16] 

according  to  the  definition  of  fx  given  by  equation  [14].  We  also  assume  that 
there  is  no  change  in  the  partial  pressure  of  solvent  vapor  during  the 
experiment,  which  is  to  say  that  the  solvent  evaporates  as  quickly  as  it  leaks 
from  the  vessel ,  as  would  certainly  be  the  case  except  when  the  leak  is 
disastrously  large  or  the  solvent  quite  nonvolatile. 

R  vessel:  We  now  consider  the  "R"  or  "reference"  vessel,  which  contains  only 
a  gas  phase.  Since  it  is  a  rigid  vessel,  the  following  pressure -temperature 
relation  is  true  for  its  contents : 

- - -  -  -  [17] 

•  _  0+  _  0+  >T>0+ 

Zr  nR  T  Zr  nR  T 

We  shall  allow  for  leakage  from  vessel  R  as  well  as  from  the  other,  and  we 
define  a  leak  coefficient  fR  as  follows: 
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[18] 


Substitution  of  equation  [18]  into  equation  [17]  and  rearrangement  of  the 
result  yields  an  expression  for  T®'*’: 


^0+ 


Zr'  ^R  T“ 


[19] 


Both  vessels:  Finally,  the  pressure  difference  between  the  two  flasks  is  known 
at  time  infinity  and  can  be  used  to  eliminate  from  the  equations: 


AP-  -  pR-  -  Px“  [20] 

Further,  at  time  0+  both  vessels  have  the  same  temperature  and  pressure 
P®'*',  and  therefore  the  same  compression  factor  2°*,  and  so  the  subscript 
denoting  which  vessel  we  are  considering  can  be  omitted.  We  now  can  obtain 
the  desired  expression  for  the  Henry's  coefficient  H  by  combining  equations 
[11],  [13],  [15],  [16],  [19],  and  [20],  which  is  straightforward: 


H" 


_ ^ 

(fx  /  H°*  ) 


(  Pr-  -  AP-  -  Pi-  -  fxPoth°~  ) 

«i  (  P®-  -  '  Potb“'  )  +  ag 


[21] 


where  -  Pl  Vl  /  Mj 
«2  -  Vg  /  RT- 

*3  “  /^r)  Pr’/Zr* 

-  (  Pr-  -  AP-  -  Pl*  -  fxP„,h®’  )/zx"  -  fx  (  Pi°‘  +  Poth°’)/z°" 


Thus,  when  the  system  is  not  ideal,  meaning  that  the  "effective"  value  of  H 
depends  on  pressure,  H.jj"  can  be  computed  using  equation  [21]  providing  that 
we  know  the  value  of  the  Henry's  coefficient  at  pressure  P°  . 

If  H  is  independent  of  pressure,  that  is,  if  the  system  is  within  the 
Henry's  law  region  at  both  pressure  P®'  and  pressure  P*  then  we  may  substitute 
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the  relation  H“  -  H°'  into  equation  [21]  and  rearrange  to  yield 


H 


«2  «3 


[22] 


where  -  (  Pr"  -  AP"  -  Pi*  )  -  f*  (  P°'  -  Pi“'  ) 
and  a2,  and  a3  are  as  defined  for  equation  [21], 

It  is  interesting  to  note  that  has  disappeared  from  the  equations 
altogether.  The  reason  is  that,  since  both  sides  of  the  apparatus  are 
elevated  quickly  to  the  same  pressure  at  time  -  0,  the  initial  differential 
pressure  is  zero.  If  the  apparatus  consisted  of  only  one  vessel,  then 
knowledge  of  the  total  pressure  at  time  -  0+  would  be  required  for  solving  for 
H.  As  it  is,  the  final  differential  pressure  (AP“)  describing  the  change  in 
pressure  of  the  system  due  to  absorption,  the  final  total  pressure  (Pr"),  and 
the  initial  total  pressure  at  time  0-  (P°‘)  are  sufficient  for  defining  the 
system. 
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III.  Computation  of  the  Ostwald  coefficient 

Henry's  law  is  commonly  encountered  in  studies  of  the  thermodynamics  of 
gas/liquid  systems.  Therefore,  equations  for  computation  of  Henry's 
coefficient  were  derived  in  the  previous  section.  However,  computation  of  the 
Henry's  coefficient  requires  knowledge  of  the  mole  fraction  of  solute  in  the 
condensed  phase,  which  is  possible  only  if  one  knows  the  average  molecular 
weight  of  the  solvent,  and  one  cannot  expect  to  determine  this  for  a  tissue 
sample.  For  tissue  samples  we  instead  compute  the  Ostwald  coefficient  F, 
defined  as  follows: 

Fj  1  -  (volume  of  dissolved  gas  "2",  evaluated  at  1  ATA  and  the 

experimental  temperature)  per  (volume  of  condensed  phase  "1"). 
Note  that  Fj  i  is  never  pressure -independent;  in  the  Henry's  Law  region  F2  ^  is 
directly  proportional  to  Pj. 

From  the  above  discussion  we  deduce  the  following  expressions  for  the 
number  of  moles  of  dissolved  gas  at  time  0-  and  time  ■»  in  our  experiment: 


(1  ATA)  Vl 

0-  „  _ _  jtO- 


R  T*  z 


1  ATA 


[23a] 


(1  ATA)  V, 

n,-  -  -  F-  [23b] 

R  T-  ATA 


Given  that  Henry's  Law  is  valid  (Fj  j  is  directly  proportional  to  P2)  at 
the  pressure  P®",  we  can  write  the  Ostwald  coefficient  at  time  0-  as 


r°'  -  r,  ata  P2°‘  /  (1  ata) 

-  Ttata  (  -  Pi®'  -  Poth°*)  /  (1  ata)  [24] 


Combining  equations  [23]  and  [24], 


n,°'  - 


ATA 


(  P®-  -  Pr  -  P^th"') 


0- 


o-> 


R  T’  z 


1  ATA 


[25] 
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If  we  assume  that  Henry's  Law  is  valid  also  at  the  final  system  pressure 
then  we  can  write  an  expression  for  the  Ostwald  coefficient  at  pressure 
Px'  that  is  analogous  to  equation  [25]; 


V,  r 


1  ATA 


n, 


R  T*  z 


(  P-  -  Pi*  -  P„th‘) 


[26] 


1  ATA 


Now,  we  can  derive  the  expression  for  the  Ostwald  coefficient  in  terms  of  the 
measured  parameters  and  variables  in  our  experiment.  The  procedure  is 
identical  to  that  used  in  section  II  except  that  nL°'  and  n^"  are  represented 
in  the  mass  balances  in  terms  of  F,  as  shown  above.  The  algebra  will  not  be 
shown  this  time. 

If  Henry's  Law  is  considered  to  apply  at  all  times  (ideal  behavior),  then 
we  can  use  equation  [26]  and  [25]  and  solve  for  F^ 


1  ATA 


Vq  «3 

Vl  «* 


ATA 


[27] 


where  Cj  and  are  defined  as  they  were  in  equation  [22]  .  The  Ostwald 
coefficient  at  any  pressure  other  than  1  ATA  is  calculated  using  the  simple 
proportionality  between  F  and  pressure  that  applies  in  the  Henry's  Law  region. 

If  there  is  non- ideal  solution  behavior  at  pressure  P^”  then  we  must  use 
equation  [23b]  rather  than  equation  [26]  to  represent  nL*,  although  we  still 
assume  the  validity  of  Henry's  Law  at  pressure  P°'.  Once  this  is  done,  we 
find  the  following  result  for  the  Ostwald  coefficient  at  Pjj“: 


F“  - 


ATA 


(  P' 


lO- 


0- 


0- 


^oth 


)  +  ATA  *3 


(1  ATA)  Vl 


[28] 
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IV.  Uncertainty  in  the  computed  gas  solubility:  Propagation  of  error 

The  partial  derivative  of  H  with  respect  to  each  par5uneter  has  been 
calculated  so  that  one  can  determine  the  sensitivity  of  the  calculated  value 
of  H  toward  the  measurement  error  in  each  parameter.  We  assume  in  the 
following  discussion  that  H  is  Independent  of  pressure,  and  so  the  derivatives 
are  evaluated  from  equation  [22]. 

Although  the  following  discussion  is  concerned  with  errors  in  the 
calculated  values  of  the  Henry  coefficient  H,  it  easy  to  show  that  the 
relative  errors  in  the  calculated  values  of  the  Ostwald  coefficient  F  are  of 
equal  magnitude,  but  have  the  opposite  sign.  In  other  words, 

1/T  (ar/dB)  -  -1/H  (dH/dB) 

where  B  is  a  parameter  subject  to  measurement  error. 


dH/dVL  -  H/Vl  [29] 

dH/dVs  -  -H/Vg  [30] 

dH/dT*  -  H/T*  [31] 

dH/dP°'  -  - ^  fx  [32] 

*2  «3 

«i 

dH/dPR-  -  -  [  «3  +  (  1/Zx’ - )  «*  ]  [33] 

*2  *3^ 

dH/dPiO-  -  - -  fx  (  «3  +  )  [34] 

*2  «3 

-  «1 

dH/dPi*  -  - -  (  «3  +  «*/2x'  )  [35] 

«2  ®3 

-  «1 

dH/dAP"  -  -  (  eg  +  «*/zx*  )  [36] 
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aH/afji 


[37] 


«2  «3^ 

«i  fx 

®2  *3^  2r* 

-  “l 

aH/afx  -  -  {  a^(9°-  -  Pi°-) 

«2  «3" 

+  [  PrV(W)  -  (Pi“'  +  Poth°‘)/2°"  +  Poth°Vzx‘  ]  }  [38] 


Note  that  our  error  analysis  is  designed  only  to  estimate  the  error  in  the 
calculated  H  resulting  from  errors  in  measuring  the  individual  parameters . 

A  truly  comprehensive  error  analysis  would  have  to  address  the  interactions 
between  parameters.  For  example,  our  representation  of  aH/ST”  is  actually 
quite  simplistic.  The  parameter  T“  appears  in  equation  [22]  only  as  a 
determinant  of  gas  pressure  according  to  an  equation  of  state,  and  therefore 
the  above  expression  for  dH/dV  accounts  only  for  the  errors  in  computing 
amounts  of  gas  that  result  from  an  error  in  T".  In  reality,  an  error  in  T" 
means  that  the  water  bath  temperature  is  different  from  its  set  point,  meaning 
in  turn  that  the  Henry's  coefficient  and  the  vapor  pressure  of  water  are 
different  from  the5r  values  at  the  set  point  temperature.  The  total  error  in 
H  resulting  from  an  error  in  T*  is  seen  to  be  a  relatively  complicated  matter 
to  analyze,  and  we  have  not  attempted  to  quantify  it.  This  error  is  minimized 
by  controlling  the  experimental  conditions  as  closely  as  possible:  the  water 
bath  temperature  is  maintained  to  within  ±  O.OS’C  of  its  set  point  during  the 
experimental  period. 

Also  note  that  equation  [37]  and  [38]  allow  us  to  calculate  errors  in  H 
due  to  the  measurement  error  of  a  leak  on  either  side  of  the  apparatus.  As 
will  be  shown  later,  actual  errors  in  H  may  be  more  sensitive  toward  the  "net 
leakage"  represented  by  (f^  -  f*)  instead. 
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Let  us  assign  plausible  values  to  the  parameters  and  get  numerical  values 
for  the  derivatives  listed  as  equations  [29] -[38].  We  want  to  see  whether  the 
experimental  results  should  be  more  precise  at  higher  or  lower  total  pressure. 
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Given:  H  -  7.6-10^  mmHg/mol  fraction  (approximately  the  average  reported 
value  for  nitrogen  in  water  at  37*0) 

Vl  -  120  cm^ 

Vg  -  30  cm^ 

T“  -  310.15*K  (37*0) 

P°'  -  760  mmHg 

Pi°*  -  Pj*  -  Pi’'“P  “  47  mmHg  (the  vapor  pressure  of  water  at  37*C) 

Poth“'  “  W  -  159.6  mmHg  (21%  of  P°-) 

“  I'O  no  leaks) 

The  results  are  expressed  as  per  cent  change  in  the  calculated  H  per  unit 
of  error  in  the  parameter  value: 

Pr*  (mmHg) 


1000 

4000 

7600 

AP* 

12.8 

173 

366 

mmHg 

an/avL 

0.833 

0.833 

0.833 

%/cm® 

an/avg 

-3.33 

-3.33 

-3.33 

%/cm® 

aH/aT“ 

0.322 

0.322 

0.322 

%/*C 

aH/ap°' 

-0.440 

-O.U326 

0.0154 

%/mmHg 

an/aPR- 

0.440 

0.0326 

0.0154 

%/mmHg 

aH/aPi°' 

8.23 

0.610 

0.289 

%/mmHg 

an/aPi* 

-8.23 

-0.610 

-0.289 

%/mmHg 

an/aAP" 

-8.23 

-0.610 

-0.289 

%/mmHg 

aH/afR 

77.9 

23.1 

20.8 

%/% 

an/afx 

-77.4 

-23.0 

-20.8 

%/% 

Clearly,  we  can  expect  higher  precision  at  higher  pressures. 

A  question  of  practical  importance  is,  how  much  water  or  tissue  should  be 
loaded  into  the  "X"  vessel  so  that  experimental  precision  is  maximized?  We 
therefore  calculate  the  partial  derivatives  for  various  amounts  of  water. 


21 


assuming  a  vessel  volume  of  150  cm®  (note  that  H  is  calculated  from  the  ratio 
of  Vl  to  Vq,  not  from  the  individual  values  of  and  Vg) . 

Given:  H  -  7.6-10^  mmHg/mol  fraction 
T“  -  310.15*K  (37'C) 

P°'  -  760  mmHg 
Pg"  -  1900  mmHg 

Pj®'  -  Pi“  -  Pi''*P  -  47  mmHg  (the  vapor  pressure  of  water  at  37*C) 

Poth°’  -  Poth“  -  159.6  mmHg  (21%  of  P®') 

“  fx  “  1-0  (i-®-.  no  leaks) 

We  calculate  the  following: 


Vl 

Vg 

-  75  cm® 
“75  cm® 

Vl  -  120  cm® 
Vg  -  30  cm® 

'  Vl  -  135  cm® 

'  Vg  “  15  cm® 

AP” 

15.9 

61.0 

129 

mmHg 

an/avL 

1.33 

0.833 

0.741 

%/cm® 

an/avg 

-1.33 

-3.33 

-6.67 

%/cm® 

an/aT” 

0.322 

0.322 

0.322 

%/*C 

an/ap®’ 

-0.0890 

-0.0927 

-0.0989 

%/mmHg 

an/aPj,” 

0.0890 

0.0927 

0.0989 

%/mmHg 

aH/aPi°* 

6.39 

1.73 

0.876 

%/mmHg 

an/aPi” 

-6.39 

-1.73 

-0.876 

%/mmHg 

an/aAP” 

-6.39 

-1.73 

-0.876 

%/mmHg 

an/afR 

120.0 

31.2 

14.8 

%/% 

an/afx  - 

117.0 

-31.0 

-15.1 

%/% 

It  appears  that  errors  from  most  sources  can  be  reduced  simply  by  using  a 
higher  liquid- to-gas  ratio.  (The  notable  exception  is  the  error  in  H  arising 
from  an  error  in  Vg;  Vg  is  calculated  by  subtracting  from  the  total  volume 
of  vessel  X.)  The  error  in  H  arising  from  a  given  absolute  error  in  AP" 
decreases  at  higher  liquid/gas  ratios  because  AP”  is  roughly  proportional  to 
Vi/Vg,  so  that  the  relative  error  in  AP”  decreases  with  increasing  V^/Vg. 
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Similarly,  3H/3Pi°‘  and  dH/dPj”  decrease  with  increasing  liquid/ gas  ratio 
because  the  effect  of  erroneous  values  of  P^®'  and  P^"  is  to  confound  the 
measurement  of  AP". 

The  liquid  volume  probably  should  not  be  set  above  120  cm^  in  the  current 
vessel,  because  the  use  of  a  greater  volume  would  place  some  of  the  liquid  in 
the  neck  of  the  vessel  (see  Figure  1).  Since  the  liquid  mass  is  agitated  only 
by  a  stirring  bar  in  the  bottom  of  the  vessel,  the  portion  within  the  neck 
could  be  poorly  mixed,  so  that  passive  diffusion  might  be  the  primary 
mechanism  of  solute  transport  through  that  portion.  This  would  necessitate  a 
longer  equilibration  time,  which  does  not  seem  desirable  since  accumulated  gas 
leakage  probably  is  roughly  proportional  to  the  duration  of  the  experiment 
(i.e.,  f]i  and  f^^  increase  with  time). 

The  numerical  values  of  the  partial  derivatives  can  be  converted  to 
errors  in  H  based  on  the  manufacturers'  claimed  precision  of  the  instruments 
we  use.  (Of  course,  the  actual  precision  in  our  experiment  could  be 
substantially  better  than  guaranteed  by  the  manufacturers  since  we  rely  on 
local  calibrations  over  relatively  narrow  ranges  of  the  measured  variables, 
whereas  the  specifications  normally  apply  to  the  full  scale.  Our  precision 
likewise  could  fail  to  meet  claimed  specifications  because  of  poor  quality 
control  or  damage  to  the  instruments.)  The  manufacturers'  specifications  for 
our  instruments  are  given  in  Appendix  A.  The  overall  estimated  error  is  the 
sum  of  the  individual  error  estimates.  We  show  this  calculation  for  one 
experimental  case : 
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Given:  H  -  7. 6' 10^  mniHg/mol  fraction 
T  -  310.15*K  (37’C) 

P“‘  -  760  mmHg 


Pr”  -  1900  mmHg 

-  Pi*  -  Pi''“*’  -  47  mmHg  (the  vapor  pressure  of  water  at  37*C) 
Poth®’  -  Poth"  -  159.6  mmHg  (21%  of  P®*) 


”  1 • ®  ( 1 • ® • f  no  leaks ) 
Vl  -  120  cm® 

Vg  -  30  cm® 


We  find  that  AP"  -  61.0  mmHg. 
error  source 


Vl  (0.833  %/cm®) 

Vg  (-3.33  %/cm®)  X 
T“  (0.322  %  /‘O  X 


X  (±  0.1  cm®) 


(±  0.1  cm®) 
(±  0.1*0 


±  0.083  % 
±  0.333  % 
t  0.032  % 


(-0.0927  %/mmHg) 
(0.0927  %/mmHg) 


X 

X 


(±  1  mmHg) 
(±9  mmHg) 


±  0.093  % 
±  0.834  % 


Pi®'  (1.73  %/mmHg) 

X 

(+  0.2  mmHg)  - 

+  0.346  % 

Pi" (-1.73  %/mmHg) 

X 

(+  0.2  mmHg)  - 

-  0.346  % 

AP*(-1.73  %/mmHg) 

X 

(±  0.031  mmHg) 

-  ±0.054 

fR  (31.2  %/%)  X  ( 

+ 

0.1  %  )  -  + 

3.116  % 

fx(-31.0  %/%)  X  ( 

+ 

0.1  %  ) 

3.105  % 

total 

+  4.89  %,  -  4.88 

% 

The  first  total  error  listed  above  was  obtained  by  adding  together  all  of 
the  possible  positive  contributions,  and  the  second  total  is  the  sum  of  all 
the  possible  negative  contributions.  One  point  to  note  is  that  leakage  and 
errors  in  the  assumed  partial  pressures  of  water  vapor  are  biased  errors. 

This  is  because  1)  gas  leaks  out  of  the  system  but  not  in,  and  2)  the  gas 
phase  in  vessel  X  can  be  below  the  dew  point  but  not  above  it.  The  remaining 
errors  presumably  are  unbiased. 

The  major  observation  from  the  above  examples  is  that  leakage  seems  to 
have  the  greatest  potential  as  an  error  source.  Based  on  preliminary 


24 


measurements  of  leakage  rates ,  the  leakage  during  an  experiment  is  taken  to  be 
0.1%  in  the  above  example,  resulting  in  an  error  in  H  of  3%  if  only  one  vessel 
leaks.  Errors  of  10%  or  more  would  result  from  a  leakage  of  only  0.5%  of  the 
contents  of  one  vessel  or  the  other,  which  is  a  very  plausible  occurrence  if 
an  0-ring  or  valve  does  not  seat  perfectly.  Note,  however,  that  equal  leakage 
from  both  vessels  produces  a  net  error  of  close  to  zero.  Obviously,  the 
difference  in  leakage  between  the  two  vessels  is  far  more  important  than  the 
total  leakage  from  both  vessels. 

The  overall  error  estimated  above  applies  to  a  single  measurement.  Our 
actual  estimate  of  precision  will  be  obtained  through  statistical  analysis  of 
our  entire  body  of  data  for  each  gas.  If  Henry's  Law  is  valid  for  a  gas  over 
the  entire  pressure  range  investigated,  then  H  equals  the  slope  of  a  plot  of 
y2^  versus  X2.  In  that  case,  we  will  use  linear  regression  to  calculate  the 
slope  and  use  the  t  distribution  to  estimate  a  confidence  interval  on  the 
slope.  If  the  F  test  indicates  statistical  support  for  modelling  the  system 
as  non-ideal,  then  we  will  use  a  non-linear  least-squares  minimization 
technique  to  fit  the  data  with  the  expression  given  at  the  end  of  part  I  of 
this  section.  The  confidence  Interval  we  obtain  statistically  will  almost 
certainly  be  more  narrow  than  the  uncertainty  interval  computed  above,  barring 
grossly  optimistic  specifications  by  the  manufacturers. 
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V.  Measurements  of  gas  leakage 

The  numerical  examples  of  the  previous  section  demonstrate  that  the 
experimental  error  attributable  to  leakage  is  determined  by  the  ”net  leakage" 
between  flasks  R  and  X.  "Net  leakage"  is  defined  as  (fj^  -  fj^)  ;  fj^  and  f^^  are 
expressed  as  fractions  and  are  defined  by  equations  [14]  and  [18].  An 
examination  of  those  numerical  results  shows  that,  to  a  close  approximation, 


dH/dfR  +  dH/fx  -  H  (fx  -  fn)  PrVAP* 


[39] 


This  can  be  proven  formally  under  the  following  assumptions  (although  the 
proof  will  not  be  shown  here) : 


fx  -  fR  -  1.0  ;  Zr-  -  zx"  -  -  1.0 


Poth  “  Poth  • 


We  can  make  the  reasonable  assumption  that  the  measurement  error  Af  is  the 
same  for  both  fR  and  fx,  i.e. 


Af  -  AfR  -  Afx  [40] 

Equation  [39]  can  be  rewritten  using  equation  [40]  to  yield  a  good  estimate  of 
the  relative  error  in  the  measured  value  of  H  as  a  function  of  the  leakages 
from  both  vessels: 


AH/H  -  Af  (fx  -  fR)  pR-  /  AP-  [41] 

Accordingly,  we  measure  (fx  -  fx)  by  pressurizing  both  vessels  simultaneously 
to  the  same  pressure  (without  a  condensed  phase  in  the  X  flask) ,  closing  the 
tubing  between  them  (by  closing  valve  V3  in  Figure  1) ,  allowing  some  time  for 
leakage  to  occur,  and  then  measuring  a)  the  pressure  difference  between  the 
vessels  (using  the  differential  transducer  T1  in  Figure  1)  and  b)  the  pressure 
in  each  vessel  (using  transducers  T2  and  T3).  To  evaluate  the  data,  we  note 
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that  if  fg  and  are  close  to  1.0  then  the  following  are  true  to  a  close 
approximation: 


(fx  -  fa) 

M  ^  pO'*' 

142] 

^R  “ 

[43] 

fx  - 

p^flnal  ^  p0+ 

[44] 

where  -  pressure  in  the  R  flask  at  the  end  of  the  experiment  (mmHg) ; 

-  pressure  in  the  X  flask  at  the  end  of  the  experiment  (mmHg) ; 

P°'*'  -  initial  system  pressure  (mmHg) ; 

^pfinai  _  pressure  difference  at  the  end  of  the  experiment, 

(P^finai  .  (nmiHg). 

Gas  solubility  measurements  can  begin  when  gas  leakage  is  deemed  acceptably 
small  based  on  the  relative  error  in  H.  An  example  calculation  is  presented; 
Given:  H  -  7,6-10^  mmHg/mol  fraction 


T"  -  310.15*K  (37*C) 

P°*  -  760  mmHg 
Pr*  -  1900  mmHg 

Pj®*  -  -  Pi'^***  -  47  mmHg  (the  vapor  pressure  of  water  at  37*C) 

Poth°'  -  Poth”  -  159.6  mmHg  (21%  of  P®*) 

"  I'O  (i-e-,  no  leaks) 

Vl  -  120  cm® 

Vq  -  30  cm® 

We  find  that  AP*  -  61.0  mmHg 

dH/dfx  -  31.2  %/% 
dH/dfR  -  31.8  %/% 

dH/dfR  +  dH/fx  -  31.1  H  (fx  •  fg)  %/% 


It  already  has  been  computed  (near  the  end  of  part  IV)  that  the  other 
error  sources  may  contribute  up  to  2.12%  to  the  relative  error  in  H  for  this 
example.  If  we  want  no  more  than  10%  total  relative  error  in  H,  then  we  will 
not  accept  any  leak  rate  that  itself  produces  more  than  7.88%  error  in  the 
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measured  H.  Accordingly,  if  it  were  known  somehow  that  only  one  of  the 
vessels  could  leak,  then  the  maximum  acceptable  leak  from  the  X  vessel  would 
be  4.8  mmHg  (or  f^  -  0.9974),  and  the  maximum  acceptable  leak  on  the  R  vessel 
would  be  4.7  mmHg  (fj^  -  0.9975).  Based  on  preliminary  leakage  measurements, 
we  can  take  the  typical  value  of  Af  (i.e.,  the  error  in  the  assumed  value  of 
f)  to  be  0.001  -  0.1%;  that  is,  when  computing  H  we  normally  would  assume  zero 
leakage,  or  f  -  1.0,  and  actual  values  of  f  are  typically  0.999.  Then,  the 
maximum  acceptable  value  of  the  net  leak  parameter  (f^  -  f^)  for  this  example 
is  found  to  be  0.00253  -  0.253%,  or  the  maximum  acceptable  differential  leak 
rate  between  the  two  vessels  is  4.8  mmHg/experimental  period. 
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Appendix  A:  Listing  of  components  In  the  experimental  apparatus 

The  following  types  of  equipment  will  be  used  in  measuring  gas  solubility: 

*  Electronic  pressure  transducers 

*  Two  glassware  flask  assemblies 

*  Thermometer  with  400 -series  cables 

*  Two  water  Baths 

*  Primary  pressure  standard  (dead  weight  tester) 

The  specific  components  obtained  for  this  project  are  as  follows: 

1.  Electronic  pressure  transducers 

A.  Differential  transducer  -  MKS,  Valley,  MD 

a.  Sensor  head: 

•  Series  170,  Type  310 

•  Pressure  range:  1000  Torr 

•  Accuracy:  ±  0.05  %  of  readings 

•  Dimensions:  2  1/2"  x  3"  x  1  3/4" 

■  Modified  In-house  for  underwater  use 

b.  Electronic  unit  -  MKS,  Valley,  MD 

•  Type  170M-6C 

•  Operation  Temperature:  15-40*0 

•  Input:  105-130  VAC 

•  Outpu'  .  0-10  VDC 

•  Dimensions:  19"  x  4  1/2"  x  7  1/4" 

B.  Absolute  Pressure  Gauges  -  Valldyne  Engineering  Corp., 

North^'ldge,  CA 

•  Type:  DP15,  diaphragm  replacement 

•  Line  pressure:  3000  pslg 

•  Pressure  Range:  variable,  depends  on  which  diaphragm 

is  used 

•  Accuracy:  0.5  %  of  full  scale 

2.  Flask  assemblies  •  MMRI  Biomedical  Engineering  Division,  Bethesda,  MD 

•  Pyrex  flask,  100  cm^ 

■  Stainless  steel  top 

3.  Thermometer  -  Yellow  Springs  Instrximents  Co.,  Yellow  Spring,  OH 

•  Model  49TA 
'  Series  400 

•  Temperature  range:  20-43*C 

•  Accuracy:  t  0.05*C  of  readings  from  30-40*C 

with  ambient  temperature  at  15-35*C 

•  Dimensions:  3/4"  x  9  1/4"  x  5  3/4" 

4 .  Water  baths 

A.  Accumulated  bath  -  NMRI  Biomedical  Engineering  Division, 

Bethesda,  MD 

’  Plexiglas 

•  Dimensions:  17"  x  13  1/2"  x  10  1/2" 


31 


B,  Supply  bath  -  Forma  Scientific  Inc.,  Marietta,  OH 

•  Model  2095 

•  Temperature  range:  -20  to  +70*C 

•  Pump  capacity:  8  GPM 

•  Sensitivity  :  (±)  0.02*C 

•  Dimensions:  26  5/8"  x  20  3/4"  X  14  1/2" 

5.  Primary  pressure  standard  -  CEC  Instrument  Inc.,  Pasadena,  CA 

•  Type  6-201-0001 

•  Cylinder  piston  type  6-001-0002 

•  Weight  set  6-002-0002  (1.5-50  psi) 

•  Pressure  range:  0.3-500  psi 

•  Accuracy:  0.015%  of  readings 

•  Dimensions:  2  1/2  x  3  x  1  3/4" 
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Appendix  B:  Measurement  of  the  volume  of  vessel  X 

The  volvime  of  our  X  flask  will  be  determined  by  measuring  the  pressure 
change  when  a  known  volume  of  gas  is  injected  into  the  vessel.  Figure  B1  is  a 
diagram  of  this  experiment.  A  steel-barrelled  syringe  of  volume  -  124.99 
cm®  (measured  gravimetrically)  is  used  to  inject  gas  through  tubing  into  the 
vessel.  The  volume  to  be  ascertained  includes  the  volume  of  the  tubing 
connecting  the  flask  with  valves  V3  and  V5  and  transducers  T1  and  T3 .  Notice 
that  the  flask,  the  valve,  and  much  of  the  tubing  are  below  the  water  surface 
whereas  the  syringe  and  some  of  the  tubing  are  above  the  surface,  in  open  air. 
The  length- to-dlameter  ratio  of  the  tubing  is  high  enough  that  heat  transport 
is  much  faster  in  the  radial  direction  than  in  the  axial  direction.  This  fact 
allows  us  to  assume  with  confidence  that,  at  steady  state,  the  portion  of  the 
system  above  water  is  entirely  at  the  ambient  temperature  while  all  of  the 
immersed  portion  is  at  the  water  bath  temperature.  A  mass  balance  around  the 
system  can  be  written  as  follows; 

**  1 

and  -  n^  +  nj  +  n,  +  nx  [B2] 

where  n  -  number  of  moles  of  air; 

total  -  subscript  denoting  the  whole  system; 

1,  2,  3,  and  X  -  subscripts  denoting  portions  of  the  system, 
consistent  with  Figure  Bl; 

0-  -  superscript  denoting  the  moment  just  before  injection; 
m  -  superscript  denoting  time  infinity. 

Ue  also  write  the  following  to  indicate  the  difference  in  system  pressure 
between  time  -  0-  and  time  -  •: 

pfinal  _  pO-  +  |g3j 
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We  measure  using  both  the  differential  pressure  transducer  T1  (see 

Figure  1)  and  pressure  transducer  T3,  both  of  which  are  referenced  to  the 
atmosphere  during  this  experiment.  Since  we  see  no  strong  reason  to  prefer 
one  of  these  transducers  to  the  other,  we  pool  the  data  from  them.  P°'  is 
equal  to  the  ambient  pressure  P.„h  because  valve  V5  is  open  to  the  atmosphere 
until  immediately  prior  to  injecting  the  air. 

Ve  denote  the  temperature  of  those  portions  of  the  system  above  water  as 
T.„k  at  time  -  0-  and  time  -  <».  The  immersed  portions  have  temperature 
at  time  -  0-  and  time  -  •».  At  intermediate  times  the  temperatures  are  higher 
because  of  the  heat  evolved  upon  compression.  Finally,  we  use  the  ideal  gas 
law  to  relate  pressure  and  temperature  at  both  time  points,  and  we  assume  zero 
leakage  from  the  system.  Combining  equations  [B1]-[B3],  we  obtain 


Vv  - 


^omb  ^bath 


‘'bath 


AP 


final 


^sflib 


where  P^  -  ambient  pressure  (mmHg) .  The  ambient  pressure  is  measured  using 
a  mercury  manometer  and  the  ambient  temperature  is  measured  using  a 
high-precision  thermometer  certified  by  the  National  Bureau  of  Standards.  We 
measure  the  volume  (V2  +  V3)  by  filling  with  water  from  a  graduated  syringe 
after  closing  the  valve  shown  in  Figure  Bl.  We  then  calculate  V2  from  the 
known  diameters  of  the  tubing  comprising  V2,  and  we  obtain  V3  by  subtracting 
V2  from  (V2  +  V3) . 
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Appendix  C:  Nomenclature  for  sections  II  and  Ill  of  "Data  analysis" 
p  -  density  (g/cm®) ; 

H  -  Henry's  coefficient  (mmHg/mol  fraction); 

M  -  molecular  weight  (g/gmol) ; 
n  -  molar  quantity  of  the  solute  (gmol); 

P  -  absolute  total  pressure  (mmHg) ; 

P  -  absolute  partial  pressure  (mmHg) ; 

AP  -  differential  pressure  (mmHg) ; 

R  -  gas  law  constant  (mmHg-cmVgniol-*K) ; 

T  -  absolute  temperature  (K) 

V  -  volume  (cm®) ; 

z  -  compressibility  factor  (dimensionless). 

subscripts  and  superscripts: 

G  -  a  subscript  that  represents  the  gas  phase; 

L  -  a  subscript  that  represents  the  liquid  phase; 

1  -  a  subscript  that  represents  the  solvent; 

2  -  a  subscript  that  represents  the  solute; 

oth  -  a  subscript  that  represents  gases  other  than  the  test  gas  or  the 
solvent  vapor; 

R  -  a  subscript  that  represents  the  "reference"  flask  (the  one  containing 
only  gas) ; 

X  -  a  subscript  that  represents  the  "experimental"  flask  (the  one 
containing  liquid  and  gas); 

0-  -  a  superscript  that  represents  the  time  just  before  pressurization; 

0+  -  a  superscript  that  represents  the  time  Just  after  pressurization; 

•e  -  a  superscript  that  represents  infinite  time  after  pressurization  in  a 
solubility  measurement. 
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